Tema 2 


Vectores y matrices 

Guillermo Peris Ripollés 


Objetivos 

Cuando finalice este tema, el alumno deberá ser capaz de: 

■ Definir vectores y matrices con Octave. 

■ Comprender y utilizar la notación de dos puntos para la creación de listas. 

■ Controlar el acceso y modificación de elementos y conjuntos de elementos de 
vectores y matrices. 

■ Realizar operaciones escalares con matrices. 

■ Conocer los distintos tipos de matrices predefinidas por Octave. 

Aplicación 

Cuando finalice este tema, el alumno deberá ser capaz de resolver problemas como el 
siguiente, cuya resolución se indica a lo largo del propio tema. 


Cinética de reacciones consecutivas 


Algunas reacciones transcurren mediante la formación de un intermediario, como 
en las reacciones consecutivas de primer orden A —> B —> C . Las ecuaciones que 
rigen estas reacciones, y que proporcionan la concentración de cada especie en 
función del tiempo t, son las siguientes: 


[A] - [A] 0 e- k ^ 

[B] = fci[A]o 


e -fcií _ g -k 2 t- 


[C] — [A] o { 1 + 


k 2 — k\ 
kie~ k2t - k 2 e~ klt 
k 2 — k\ 


donde [A], [B] y [C\ son las concentraciones de las respectivas especies, [A]q 
la concentración inicial de A y k\ y k 2 las constantes de velocidad. Escribe un 
programa que pida al usuario los valores de [A]o, fci y k 2 , y calcule las concentra¬ 
ciones de las tres especies en función del tiempo. Comprueba el funcionamiento 
del programa cuando [d.] 0 =1 mol/l, k-\ =0.1 s _1 y k 2 =0.2 s -1 . 
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2.1 Introducción 
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2.1. Introducción 

Hasta ahora únicamente hemos trabajado con números simples (más formalmente, 
escalares), pero ya hemos comentado que la enorme potencia de cálculo de Octave se 
encuentra fundamentalmente en el cálculo vectorial y matricial. En esta unidad se van 
a introducir estos nuevos conceptos de forma simple, aplicándolos en su parte final a 
la resolución de problemas relacionados con la ingeniería. 

2.2. Vectores 

En Octave los vectores se pueden crear introduciendo una lista de valores separados 
por espacios o comas y encerrados entre corchetes. Veamos un ejemplo a continuación: 


>t = [4 8 -2 3 5] 
t = 

4 8 -2 3 5 


En numerosas ocasiones, nos interesarán listas de valores en las que sus elemen¬ 
tos guarden una cierta estructura, relación u orden. Por ejemplo, podríamos estar 
interesados en un vector con los enteros comprendidos entre 0 y 10: 


»t = [0 123456789 10] 
t = 

0123456789 10 


Los valores de los elementos de t los hemos introducido uno a uno porque no 
son muchos pero, ¿y si hubiéramos querido introducir una lista de valores de 0 a 
100?. Para facilitar esta tarea, Octave introduce la notación de dos puntos(:). Si 
escribimos dos números enteros separados por dos puntos, Octave genera todos los 
enteros comprendidos entre ellos. Así, podríamos crear el vector t como sigue: 


»t = [0:10] ; 


Es decir, la orden [i : j ] crea el vector [i i+1 i+2 ... j-2 j-1 j ]. Si quisiéra¬ 
mos que el intervalo entre los elementos fuera distinto de 1, utilizaríamos tres números 
separados por siendo el número central el incremento: 


>s = [0:2:10] 
s = 

0 2 4 6 8 10 


En este ejemplo, el vector creado contiene los valores entre 0 y 10 separados por 2 
unidades. Este valor de incremento puede poseer cualquier valor entero, real e incluso 
negativo: 


»s = [10:-2:0] 
s = 

10 8 6 4 2 0 
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Gctave también tiene definida la función linspace (inicial, final .elementos) 
que permite crear un número determinado de valores (elementos) distribuidos homo¬ 
géneamente entre dos límites (inicial y final) . Dicho de otra forma, en lugar de 
especificar el intervalo entre elementos (como se hacía con la notación de dos puntos) 
se indica el número de elementos que queremos. Por ejemplo, 


3>s = linspace(0,10,6) 
s = 

0 2 4 6 8 10 


crea un vector de 6 elementos equiespaciados entre 0 y 10. También en este caso se 
podría crear un vector de valores decrecientes cambiando el orden de los valores inicial 
y final: 


3>s = linspace(10,0,6) 
s = 

10 8 6 4 2 0 


Ejercicios - 

@ ► 1 ¿Qué vectores crean las siguientes instrucciones? Piénsalo antes de obtener el resultado 

con Octave. 


a) 

[1:9] 

e) 

[1 

3:9] 

b) 

[-1:9] 

f) 

[9 

-2:1] 

c) 

[1:2:9] 

g) 

[9 

-3:1] 

d) 

[1:0.5:9] 

h) 

[9 

-3:-1] 


lít H ► 2 Utiliza la función linspace para crear los vectores del ejercicio anterior. 


2.2.1. Acceso y modificación de elementos de vectores 

En ocasiones, nos interesa acceder a elementos determinados de un vector previamente 
definido, o bien modificar dichos elementos. Para ello, Octave permite el uso de índices 
(expresados entre paréntesis), de forma que t(l) sería el primer elemento de t, t(2) 
el segundo, etc. Veamos un ejemplo: 


»t =[6 1 -1 12 4.4 -2.36 2 0.87 48 9 0] ; 

»t(l) 

ans = 

6 

»t(6) 
ans = 

-2.36 


También puede utilizarse la notación de dos puntos para acceder a un conjunto de 
elementos, en cuyo caso el resultado es un nuevo vector. 
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>r = 

t(1:3) 


r = 

6 

1 -1 


2> s = 

t(10:-1:7) 


s = 

9 

48 0.87 

2 


En el ejemplo anterior, el vector r incluye los elementos de las posiciones 1-2-3 de 
t, y el vector s los elementos de las posiciones 10-9-8-7, en ese orden. De hecho, es 
posible acceder a elementos cualesquiera y en cualquier orden del vector mediante el 
uso de una lista de elementos. Por ejemplo, en las siguientes líneas se crea un nuevo 
vector u eligiendo los elementos de las posiciones 3, 8, 1 y 6 de t, en este orden: 


3>u = t( [3 8 

1 6]) 

u 


: 0.87 

6 -2.36 


Mediante esta selección de elementos de un vector, también podemos modificar 
vectores ya creados, mediante órdenes de asignación. A continuación se muestra un 
ejemplo a partir del vector t: 



Fijémonos en que en algunas de las instrucciones anteriores se realizan asignacio¬ 
nes múltiples. Así, en la segunda instrucción se asignan simultáneamente 3 valores a 
3 posiciones del vector t, mientras que en la tercera se realizan 5 asignaciones. Obvia¬ 
mente, si el número de elementos a ambos lados del operador de asignación es distinto, 
Octave no realizará ninguna asignación, indicándolo con un error. 


»t(4:6) = [5 4 3 2] 

error: A(I) = X: X must be a scalar or a vector with the same 
length as I 
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Ejercicios - 

H ► 3 A partir del vector v = [1 3 5 7 9 11 13 15 17 19], ejecuta las siguientes órdenes 
en Octave, pensando previamente cuál va a ser el resultado. 

a) v(3) d) v(10: -2 :1) 

b) v(3:5) e) v([l 3 5 7 9]) 

c) v(5:-l:3) f) v(v(l:5)) 

m ► 4 Piensa cuál de las siguientes asignaciones es correcta, indicando en ese caso cuál sería 
el resultado. Utiliza el vector v definido en el ejercicio anterior. Compruébalo con Octave. 

a) a = v(3) c) b(l:2) = v(7:9) 

b) b = v(3:5) d) b(l:3) = v(7:9) 


2.3. Matrices 


Denominamos matriz de dimensiones m x n a un conjunto de números ordenados en 
m filas y n columnas. 


Olí 

012 • 

Oln 

021 

022 • 

’ • 02n 

®ml 

Om2 

o mn 


Un vector no es más que un caso particular de matriz en el que una de sus dimen¬ 
siones es 1. Denominaremos vector fila a una matriz 1 x n y vector columna a una 
matriz m x 1. 


Vector fila [ai a 2 • • • a n ] 


Vector columna 


ai 

«2 


a 


m 


Para introducir matrices en Octave utilizaremos corchetes (al igual que con vecto¬ 
res), separando las filas por puntos y coma (;) y las columnas de cada fila con espacios 
(o comas). Las filas también pueden separarse con cambios de línea. A continuación se 
muestra la definición de una matriz A de dimensiones 3 x 4 de dos formas distintas 0 



1 ¿Se te ocurre una forma más rápida de introducir esta matriz? 
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»A = 

[12 3 4 

5 6 7 8 


A = 

9 10 11 

12] 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


La función size aplicada sobre una matriz (o vector) devuelve un vector de dos 
componentes con sus dimensiones en filas y columnas: 


3> size(A) 


ans= 

3 4 


Podemos acceder a los elementos individuales de esta matriz indicando entre pa¬ 
réntesis la fila y columna del elemento separadas por comas. Así, 


»A(2,3) 
ans = 

7 


De nuevo, la potente notación de dos puntos de Octave puede utilizarse para crear 
subnratrices a partir de una matriz original. Considera los siguientes ejemplos, en los 
que se utiliza la matriz A definida anteriormente: 


»b = 

b = 

A(2:3,3) 


7 



11 



»c = 

A(3,:) 


c = 



9 

10 11 

12 

>d = 

d = 

A ( [2 3] , [2 

4]) 

6 

8 


10 

12 



En la primera expresión, se crea un nuevo vector columna b tomando de las filas 2 
y 3 los elementos de la tercera columna. En la segunda orden, utilizamos el operador 
sin números que lo delimiten para indicar que nos interesan todos los elementos de 
la fila 3 de A. Por último, en la expresión final creamos una nueva matriz tomando los 
elementos de las posiciones 2 y 4 de las filas 2 y 3. 

Al igual que se vió en el caso de vectores, también podemos modificar los elementos 
de una matriz con la instrucción de asignación =, siempre que las dimensiones a ambos 
lados del signo igual sean las mismas: 
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>A([2 

3] , 

[2 4]) = [0 1; 2 3] 

A = 

1 2 

3 

4 

5 

0 

7 

1 

9 

2 

11 

3 


Por último, también podemos construir nuevas matrices yuxtaponiendo matrices 
ya existentes, siempre que éstas tengan las dimensiones adecuadas. Fíjate en los si¬ 
guientes ejemplos: 


>X = 

[i 

2 ; 3 4] 

X = 

1 

2 


3 

4 


»Y = 

[5 

6 ; 7 8] 

Y = 

5 

6 


7 

8 


»Z = 

[3 

4 5 6 ; 7 8 9 10] 

z = 

3 

4 

5 6 

7 

8 

9 10 

»B = 

[X 

Y] 

B = 

1 

2 

5 6 

3 

4 

7 8 

»C = 

[X 

; Y] 

c = 

1 

2 


3 

4 


5 

6 


7 

8 


»D = 

[X 

Z] 

D = 

1 

2 

3 4 5 6 

3 

4 

7 8 9 10 

»E = 

[X 

; Z] 

error : 

number of columns must match (4 != 2) 


Ejercicios - 

S ► 5 Define con Octave las siguientes matrices y comprueba sus dimensiones con la función 

size. 


A = 


1 

7 


2 3 
8 9 


4 

10 





'1 

7' 





2 

8 


5 

6' 

B = 

3 

9 

C = 

11 

12 

4 

10 




5 

11 





6 

12 



0.6 

1.5 

2.3 

-0.5 

8.2 

0.5 

-0.1 

-2.0 

5.7 

8.2 

9.0 

1.5 

0.5 

0.5 

2.4 

0.5 

1.2 

-2.3 

-4.5 

0.5 
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@ ► 6 ¿Cuál será el resultado de las siguientes expresiones en Octave? Comprueba tus afir¬ 
maciones con el propio programa. Utiliza las matrices definidas en el ejercicio anterior. 


al = A(2,3) 

i) c4 = C(l:2:5,:) 

a2 =A(3,2) 

j) c 5 =C([5 2 1], 3:-l:2) 

bl = B(2,3) 

k) D = [4:9; 1:6] 

b2 =B(3,2) 

1) E = [D A] 

el = C(2,3) 

m) A(:,2:3) = B(l:2,:) 

c2 =C(3,2) 

n) A(l,4:6) = C(4,l:3) 

a3 = A(:,2) 

ñ) B(4:6,2) = B(4:6,l) 

c3 = C(4:5,l:3) 

o) C(l,:) = A(l,l:4) 


2.4. Matrices predefinidas 

La definición de matrices mediante la introducción elemento a elemento es costosa 
cuando se trata de matrices de gran tamaño. Para facilitar la labor de crear matrices, 
Octave dispone de una serie de funciones predefinidas: 

■ Matriz de ceros : La función zeros genera una matriz que contiene O’s en 
todos sus elementos. Si la función sólo recibe un argumento (zeros (m) ) se genera 
una matriz cuadrada de m filas, mientras que si se le pasan dos (zeros(m,n)) 
se genera una matriz rectangular de dimensiones m x n. En este último caso, 
también puede utilizarse la salida de la función size para crear una matriz de 
ceros cuyas dimensiones vendrán dadas por otra matriz. Todos estos conceptos 
se entenderán mejor con los siguientes ejemplos: 


>A = 

zeros 

(3) 

A = 

0 

0 

0 

0 

0 

0 

0 

0 

0 

»B = 

zeros 

(3,2) 

B = 

0 

0 


0 

0 


0 

0 


>c = 

[ 1 2 

3 ; 4 2 5] 

c = 

1 

2 

3 

4 

2 

5 

»D = 

zeros(size(C)) 

D = 

0 

0 

0 

0 

0 

0 


■ Matriz de unos : Exactamente igual que la función zeros, pero ahora se genera 
una matriz con elementos unidad. El nombre de la función es ones. 
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3>A = ones(3,2) 
A = 

1 1 
1 1 
1 1 


■ Matriz identidad : La función eye crea una matriz en la que todos los ele¬ 
mentos son ceros a excepción de los elementos diagonales (aquellos en que los 
índices de filas y columnas son iguales) que son unos. A esta matriz se le conoce 
como matriz identidad y se denota matemáticamente como I. Los argumentos 
de esta función son similares a los de las funciones anteriormente definidas. Aun¬ 
que, formalmente, las matrices identidad son cuadradas, esta función de Octave 
permite la creación de matrices identidad rectangulares. 


>A = 

eye(3) 

A = 

1 

0 

0 

0 

1 

0 

0 

0 

1 

»B = 

eye (3 

,2) 

B = 

1 

0 


0 

1 


0 

0 


>C = 

[ 1 2 

3 ; 4 2 5] 

c = 

1 

2 

3 

4 

2 

5 

>D = 

eye(size(C)) 

D = 

1 

0 

0 

0 

1 

0 


■ Matriz diagonal: Se denomina así a aquella matriz en la que todos los elemen¬ 
tos extradiagonales son ceros, es decir, únicamente presentan valores distintos 
de cero los elementos de la diagonal. Octave nos permite crear una matriz dia¬ 
gonal a partir de un vector, que sitúa en la diagonal de la matriz. Para ello 
utilizaremos la función diag como sigue: 


>X = [4 3 2 1] ; 
»Y = diag(X) 

Y = 

4 0 0 0 

0 3 0 0 

0 0 2 0 

0 0 0 1 


Por el contrario, si el argumento de la función es una matriz, diag extrae el 
vector que contiene los elementos diagonales: 
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»A = 

A = 

[ 1 

2 3; 4 5 6; 7 8 9] 

1 

2 

3 

4 

5 

6 

7 

8 

9 

2> diag(A) 

ans = 

1 

5 

9 



■ Matriz aleatoria: La función rand(n) crea una matriz de dimensión nxn con 
elementos aleatorios distribuidos uniformemente entre 0 y 1, mientras que con 
dos argumentos (rand(m,n)) crea una matriz de dimensiones m x n. 


2>A = rand(4,5) 

A — 




0.722684 

0.678103 

0.228702 

0.648577 

0.124333 

0.336018 

0.976399 

0.458428 

0.525437 

0.069060 

0.651192 

0.483592 

0.744526 

0.952700 

0.229671 

0.464437 

0.166825 

0.163426 

0.568459 

0.038776 


Ejercicios - 

H ► 7 Haciendo uso de las funciones vistas hasta ahora, define las siguientes matrices: 


'1 

0 

0 

0' 


'7 

0 

0 


0 ' 

0 

1 

0 

0 

c) 

0 

9 

0 


0 

0 

0 

1 

0 

0 

0 

11 


0 

0 

0 

0 

1 


0 

0 

0 

13 

'1 

0 

0 

0' 


'2 

1 

1 

1' 


0 

1 

0 

0 

d) 

1 

2 

1 

1 


0 

0 

1 

1 

1 

1 

2 

1 


0 

0 

1 

1 


1 

1 

1 

2 



H ► 8 Dada una matriz A cualquiera, ¿a qué es igual la operación diag(diag(A))? Para 
comprobar tu respuesta, aplícalo a una matriz A aleatoria de dimensiones 3x3 creada con la 
función rand. 


2.5. Operaciones con matrices 

2.5.1. Transposición 


La matriz transpuesta de una matriz A de dimensión m x n es una matriz de n x m 
en la que el elemento ( i,j ) ha pasado a la posición (j, i), o dicho coloquialmente, se 
han intercambiado las filas y columnas de A. Por ejemplo, la matriz transpuesta de 


12 3 4 
5 6 7 8 


es 


1 5 

2 6 

3 7 

4 8 


Ingeniería Química 


Programación en Octave 




















2-12 


Vectores y matrices 


En Octave, la transposición de una matriz se efectúa añadiendo el apostrofe (’) al 
nombre de la matriz El 


»A= [1 2 3 4; 5 6 7 8]; 
»B = A’ 

B = 

1 5 

2 6 

3 7 

4 8 


Observamos que con la transposición se convierten vectores fila en vectores colum¬ 
na y viceversa. Esta característica puede ser muy útil para mostrar vectores en forma 
de tabla. Imaginemos, por ejemplo, que de un experimento de cinética de reacción con 
medidas de conductividad obtenemos un vector t de tiempos, conc de concentracio¬ 
nes y rho de conductividades. Podemos visualizar una tabla de valores del siguiente 
modo: 


»t = [1:5] ; 

3>conc = [0.5 0.3 0. 

05 0.01 0.00005]; 

»rho = [27.6 15.4 2 

.89 0.56 0.026]; 

[t ; conc’ rho ; ] 

ans = 

1.0000 0.50000 

27.6000 

2.0000 0.30000 

15.4000 

3.0000 0.05000 

2.8900 

4.0000 0.01000 

0.5600 

5.0000 0.00005 

0.0260 


2.5.2. Aritmética escalar 

En Octave se denominan operaciones escalares a aquellas que actúan individualmente 
sobre cada uno de los elementos de una matriz o vector. Forman parte de esta categoría 
las operaciones mixtas binarias en las que intervienen una matriz y un escalar. Por 
ejemplo, si queremos multiplicar por un mismo valor todos los elementos de una 
matriz, o bien sumarles dicho valor, realizaríamos las siguientes operaciones: 


»A 

= [123 

4; 5 6 7 8]; 

>B 

= 2*A 


B = 



2 

4 6 

8 

10 

12 14 

16 


»c 

= B 

- 2 


c = 

0 

2 

4 

6 

8 

10 

12 

14 


2 En el teclado el apostrofe suele estar en la tecla a la derecha del 0. 
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La matriz B se construye multiplicando por 2 cada uno de los elementos de A, 
mientras que C se define restando 2 a los elementos de B. 

Pero, ¿qué ocurre si queremos multiplicar dos matrices de dimensiones semejantes 
elemento a elemento?. En ese caso no podemos utilizar el operador ya que éste 
realizaría un producto matricial (explicado en una sección posterior), cuyo resultado 
sería distinto al esperado. En general, para realizar operaciones escalares con matrices 
utilizaremos los operadores antecedidos de un punto. Para entenderlo mejor, vamos a 
seguir una serie de ejemplos utilizando vectores: 


>A = [2 5 

6] ; 


>B = [2 3 

5] ; 


>C= A. *B 



C = 



4 15 

30 


>D = A./B 



D = 



1.0000 

1.6667 

1.2000 


Las sumas y restas de matrices (siempre que sus dimensiones sean las mismas) 
se consideran en el algebra matricial como operaciones elemento a elemento (al con¬ 
trario que las multiplicaciones y divisiones) por lo que, aunque podemos utilizar los 
operadores . + y . el punto no es necesario. 


>A + 

B 


ans = 



4 

8 

11 

>A - 

B 


ans = 



0 

2 

1 


También pueden realizarse exponenciaciones elemento a elemento, como se muestra 
en el siguiente ejemplo: 


»C= A. a 2 


C = 


4 25 

36 

>D= 3. a A 


D = 


9 243 

729 

»E = A. a B 


E = 


4 125 

7776 


Así pues, en general, sólo será necesario anteceder el operador con un punto en 
el caso de multiplicaciones, divisiones y exponenciaciones cuando se pretenda realizar 
una operación elemento a elemento. 
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Estas operaciones, además de ser válidas para vectores, también lo son para ma¬ 
trices, como se observa en los siguientes ejemplos: 



Ejercicios - 

0 ► 9 Calcula los valores del vector C tras ejecutar las órdenes que se indican, utilizando los 
vectores A = [2 -1 5 0] y B = [3 2-14]. Comprueba tus respuestas con Octave. 

a) C = B + A -3; 

b) C = A./B; 

c) C = A. A B; 

d) C = 2*A + A. A B; 

e) C = 2. a B + A; 

f) C = 2*B/3.*A; 

¿i H ► 10 Crea un vector x con los siguientes elementos: 

a) 1, 1/2, 1/3, 1/4, 1/5, 1/6 

b) 0, 1/2, 2/3, 3/4, 4/5 


2.5.3. Aplicación de funciones sobre matrices 

Las funciones matemáticas de Octave actúan sobre las matrices del mismo modo que 
las operaciones escalares, es decir, se aplican elemento a elemento. Veamos un ejemplo: 


x = pi* [0 : 0.1: 0.9] ; 

3>y = eos (x) 

y = 

Columns 1 through 5: 

1.0000e+00 9.5106e-01 8.0902e-01 5.8779e-01 3.0902e-01 


Columns 6 through 10: 

6.1230e-17 -3.0902e-01 -5.8779e-01 -8.0902e-01 -9.5106e-01 
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3>z = abs(y) 
z = 

Columns 1 through 5: 

1.0000e+00 9.5106e-01 8.0902e-01 5.8779e-01 3.0902e-01 
Columns 6 through 10: 

6.1230e-17 3.0902e-01 5.8779e-01 8.0902e-01 9.5106e-01 


Así, el elemento i del vector y contiene el coseno del elemento i del vector x, y el 
de z su valor absoluto. 

Ejercicios - 

11 Dado el vector t = [1:0.2:2], calcula las siguientes expresiones matemáticas con 

Octave: 

a) ln(2 + t + t 2 ) d) arctg(t) 

b) e 4 (l + cos(3í)) e) cotg(t) 

c) cos 2 (t) + sin 2 (í) f) sec 2 (i) + cotg(t) — 1 


Además, Octave dispone de una serie de funciones que actúan de forma específica 
sobre vectores y matrices. Vamos a considerar aquí algunas de ellas, de interés en 
cálculos estadísticos. En particular, vamos a destacar el uso de funciones para el cálculo 
de máximos y mínimos, sumas y productos de elementos, medias y ordenaciones. La 
actuación de estas funciones depende de que sus argumentos sean vectores y funciones, 
tal y como se explica en la Tabla |2.1| 


Función 

Aplicación sobre 

Vectores 

Matrices 

max(x) 

Valor mayor en x 

Vector fila con elemento máximo 
de cada columna 

min(x) 

Valor menor en x 

Vector fila con elemento mínimo 
de cada columna 

length(x) 

Número de elementos en x 

Dimensión máxima de la matriz 

sum(x) 

Suma de los elementos de x 

Vector fila con la suma de los ele¬ 
mentos de cada columna 

prod(x) 

Producto de los elementos de x 

Vector fila con el producto de los 
elementos de cada columna 

mean(x) 

Media de los elementos de x 

Vector fila con la media de los ele¬ 
mentos de cada columna 

std(x) 

Desviación estándar de los ele¬ 
mentos de x 

Vector fila con la desviación es¬ 
tándar de los elementos de cada 
columna 

median(x) 

Mediana de los elementos de x 

Vector fila con la mediana de los 
elementos de cada columna 

sort(x) 

Ordena el vector x de forma as¬ 
cendente 

Devuelve la matriz con los ele¬ 
mentos de cada columna ordena¬ 
dos de forma ascendente 


Tabla 2.1: Algunas funciones aplicables sobre vectores y matrices. 
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Utilicemos vectores y matrices concretos para entender mejor la acción de estas 
funciones. 


»v = 

[3 

4 1 

9 

5 6] 



v = 







3 4 

1 

9 

5 

6 



»M = 

[0. 

6 1. 

5 

2.3 

-0.5 

; 8.2 0.5 -0.1 -2.0; 5.7 8.2 9.0 1.5 ; ... 

0.5 0. 

5 2 

.4 0 

.5 

; i 

.2-2. 

.3 -4.5 0.5] 

M = 







0.6 

1. 

5 

2. 

3 - 

•0.5 


8.2 

0. 

5 - 

■0. 

1 

-2 


5.7 

8. 

2 


9 

1.5 


0.5 

0. 

5 

2. 

4 

0.5 


1.2 

-2. 

3 - 

•4. 

5 

0.5 


3>max 

(v) 






ans = 

9 






3>max 

(M) 






ans = 







8.2 

8.2 

; 9 

1 

..5 




La función max devuelve, al ser aplicada a un vector, el valor máximo del mismo, 
en este caso 9. Sin embargo, al tener como argumento una matriz devuelve un vector 
fila en el que cada elemento es el máximo de la columna correspondiente de la matriz 
original. La función sum presenta un comportamiento similar con vectores y matrices, 
sumando las componentes de un vector y las componentes de cada columna en una 
matriz, como se muestra en el siguiente ejemplo: 


3> sum(v) 
ans = 28 
3> sum(M) 
ans 

16.2 8.4 9.1 0 


La función length calcula la longitud del vector, o sea, el número de elementos que 
contiene. En el caso de una matriz, devuelve el número máximo de filas o columnas, 
es decir, la dimensión de la matriz de mayor valor. 



La función mean calcula la media de los elementos del vector, es decir, la suma de 
sus elementos dividida por el número de ellos. Esto se aplica sobre las columnas al 
tener como argumento una matriz. 


Universitat Jaume I 


Guillermo Peris Ripollés 







2.5 Operaciones con matrices 


2-17 


3>mean(v) 
ans = 

4.6667 
3>mean(M) 
ans = 

3.24 1.68 1.82 0 


Por último, la función sort ordena los elementos de un vector, mientras que en el 
caso de matrices ordena de forma independiente cada una de sus columnas. 


3> sort(v) 
ans = 

1 3 4 5 6 9 

3> sort(M) 


ans = 


0.5 

-2.3 

-4.5 

-2 

0.6 

0.5 

-0.1 

-0.5 

1.2 

0.5 

2.3 

0.5 

5.7 

1.5 

2.4 

0.5 

8.2 

8.2 

9 

1.5 


Ejercicios 


12 Determina las matrices devueltas por las siguientes llamadas a función, verificando 
posteriormente tus respuestas con Octave. Deberás definir previamente las siguientes variables: 


w = [0 3 -2 7] x= [3 


-1 5 7] 



7 

4 

-2 


a) max(w) 

b) min(y) 

c) mean(y) 

d) median(w) 


e) sort(2*w+x) 

f) sort(y) 

g) length(w) 

h) length(y) 


13 En Octave, ¿es equivalente mean(x) y sum(x)/length(x)? Compruébalo con los si¬ 
guientes vector y matriz y trata de dar una explicación a tus observaciones. 

™=[0 3 —2 7] y=^ 2 l l 

¡£¡ ► 14 ¿Cómo calcularías el valor máximo de entre todos los elementos de una matriz? 
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2.6. Aritmética matriz-matriz 

Vamos a tratar ahora las operaciones internas entre matrices, que son operaciones ma¬ 
temáticas en que ambos operadores son matrices. Aunque la suma y resta de matrices 
son operaciones elemento a elemento, las repasamos aquí por encontrarse dentro de 
las operaciones entre matrices. 


2.6.1. Suma y resta 


La suma de dos matrices únicamente es posible si éstas tienen las mismas dimensiones. 
Es decir, la función size aplicada a ambas ha de proporcionar el mismo resultado. 
La matriz suma es una matriz de las mismas dimensiones en la que cada uno de sus 
elementos es igual a la suma de los elementos que ocupan la misma posición en las 
matrices sumando. Por ejemplo, 


2 1 
4 6 


+ 


4 2 
0 1 


6 3 
4 7 


En Octave, esta suma se llevaría a cabo como sigue: 


>A = 

[2 1 ; 

4 6] ; 

»B = 

[4 2 ; 

0 1] ; 

»C = 

C = 

A + B 


6 

3 


4 

7 



La resta se define y trata en Octave de forma equivalente. 


2.6.2. Producto escalar 


Podemos definir el producto escalar entre un vector fila X y un vector columna Y de 
la misma longitud como sigue: 


X Y = [xix 2 x 3 ...x n \ 


V i 
V2 
V3 


n 

= Xi • yi + x 2 • V2 + • • • + X n • y n = ^ Xj • . 

|=i 


:Vn 


Por ejemplo, para calcular el producto escalar entre los vectores A = [2 5 -2] y 
B = [4-3 2] ejecutaríamos el siguiente conjunto de órdenes: 


»A = 

[2 

5 -2]; 

>B = 

[4 

-3 2] ; 

> A*B’ 



ans = 



-11 




El resultado puede obtenerse fácilmente a mano (2*4 + 5 * (—3) + (—2) *2 = —11). 
Fijémonos en que A es un vector fila, y que hemos tenido que transponer el vector 
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antes del producto para convertirlo en un vector columna. Si no lo hubiéramos hecho, 
habríamos obtenido un mensaje de erroi0 


> A*B 

??? Error using ==>* 

Inner matrix dimensions must agree. 


El producto escalar de dos vectores es muy importante en algebra vectorial y en 
física. 

Ejercicios - 

15 A partir de los siguientes vectores, indica cuáles de las siguientes operaciones son 
correctas, calculando su valor. 




A = [5 3 -1] 

B = 

ones( 1,3) C = [9 4 0] 



D = eye(l, 4) 

E=[-1 3 -5 7] 

a) 

A*B 



e) B*C’ 

b) 

A*B ’ 



f) B*D ’ 

c) 

B*A’ 



d) 

> 

* 

O 



g) E*D ’ 


2.6.3. Multiplicación 

La multiplicación de matrices no es una operación elemento a elemento, es decir, no 
se efectúa multiplicando los elementos correspondientes de las matrices operando. Si 
realizamos la multiplicación C = A x B, el elemento Cy no se calcula como x bij, 
sino como el producto escalar de los vectores correspondientes a la fila i de la matriz A 
y la columna j de la matriz B. En la figura [2TT| se muestra gráficamente esta operación 
[^] Aquí se observa que si queremos multiplicar dos matrices A (m x n) y B (p x q), 
el número de columnas de A debe coincidir con el número de filas de B, es decir, n = p. 


C A X B 


A 


m 



m 


fila i 


x 


columna j 


A 


P 


Y 


Figura 2.1: Multiplicación de matrices. 


Una forma sencilla de decidir si podemos multiplicar dos matrices es escribir las 
tallas de las dos matrices juntas. De esta forma, si las dos dimensiones interiores coinci¬ 
den, las matrices se pueden multiplicar y la dimensión de la matriz resultante coincide 

3 Podríamos haber obtenido el resultado correcto utilizando la siguiente expresión: sum(A.*B). 

4 ¡Cuidado!. De esta definición se infiere que el producto de matrices no es conmutativo, es decir, 
Ax B A B x A. 
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con las dimensiones externas. Por ejemplo, supongamos que queremos multiplicar una 
matriz de 2 x 3 por otra de 3 x 3 (en ese orden). 

2x3,3 x3' 

Al ser los dos números internos iguales, la multiplicación es válida, y las dimen¬ 
siones de la matriz resultante son los números externos, es decir, 2x3. Consideremos 
las matrices siguientes: 


A = 




" 

2 

5 1 ' 

B = 

0 

3 -1 





1 0 2 
-1 4 -2 
5 2 1 


El producto de estas matrices se realizaría como sigue: 


C = AxB = 


2 5 1 

0 3-1 


1 0 2 
-1 4 -2 
5 2 1 

2 • 1 + 5 • (— 1 ) + 1-5 2 -0 + 5 -4 + 1 -2 2 • 2 + 5 • (-2) + 1 • 1 

0 • 1 + 3 • (— 1 ) + (— 1 ) • 5 0 • 0 + 3 • 4 + (— 1 ) • 2 0 • 2 + 3 • (-2) + (-1) • 1 

2 22 -5 
-8 10 -7 


En Octave, esta operación se realizaría de forma sencilla con el operador de multipli¬ 
cación: 


»A 

= 

[2 5 

i; 

>B 

= 

[1 0 

2; 

»c 

= 

A*B 


C = 

2 

22 

-5 


-8 

10 

-7 

>D 

= 

B*A 



0 3-1]; 

-1 4 -2; 5 2 1] 


??? Error using ==>* 

Inner matrix dimensions must agree. 


Fijémonos en que Octave no permite el producto B x A por no ser adecuadas las 
dimensiones de las matrices. 

Recordemos que si A es una matriz, A. A 2 es la operación que eleva al cuadrado cada 
elemento de la matriz. Si queremos calcular propiamente el cuadrado de la matriz, 
es decir, A*A, podemos ejecutar la orden A a 2. De este modo, A a 4 es equivalente a 
A*A*A*A. Si tenemos en cuenta que en el producto de matrices el número de columnas 
del primer operando ha de ser igual al número de filas del segundo, deducimos que 
únicamente podemos obtener potencias de matrices cuadradas. 


>B a 3 



ans = 



27 

24 

18 

-87 

24 

-66 

33 

54 

3 
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Fijémonos en que el resultado es distinto al que se obtiene con el operador punto 
correspondiente: 


>B. a 3 



ans = 



i 

0 

8 

-i 

64 

-8 

125 

8 

1 


Ejercicios 


í& H ► 16 Calcula los siguientes productos de matrices, pensando previamente cuál va a ser el 
resultado. 



2 

1 ' 


' 1 3' 

-1 5 


' 3 

2 ' 

A = 

0 

-1 

B = 

C = 

-1 

—2 


3 

0 



0 

2 


D=[ 1 2] 


a) A*B 

b) A*C 

c) A*C ’ 


d) D*B 

e) D*B a 2 

f) (C*B)*D’ 


g) B*C ’ 

h) A*D 


2.6.4. Determinantes 

El cálculo del determinante de una matriz cuadrada se realiza con la función det. 
Los determinantes son números calculados a partir de los elementos de una matriz, 
con amplias aplicaciones en cálculos matemáticos, físicos y de ingeniería. Matemáti¬ 
camente, el determinante de una matriz A se denota como |A|. 


»A = [1 0 2; -1 4 -2; 5 2 1]; 
3> det (A) 
ans = 

-36 


2.6.5. Matriz inversa 

La inversa de una matriz cuadrada A es aquella matriz que, multiplicada por la 
matriz original, proporciona la matriz identidad, es decir, A x A -1 = A~ l x A = II . 
Se puede demostrar que únicamente tienen inversa aquellas matrices cuadradas cuyo 
determinante es distinto de cero. 

Por ejemplo, consideremos la siguiente definición de dos matrices A y B: 
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1.5000 

-0.5000 

-2.0000 

1.0000 

> A*B 


ans = 


1 0 


0 1 



>B*A 

ans 

: o 

o i 


Observamos que A ■ B = B • A y que este producto es igual a la matriz identidad 
de dimensión dos, por lo que A es la matriz inversa de B, y viceversa. A continuación 
comprobamos esta afirmación con Octave , calculando la inversa de A, comprobando 
previamente que el valor de su determinante no es cero. 


3> det (A) 
ans = 

2 

>inv(A) 
ans = 

1.5 -0.5 

-2 1 


Si definimos una matriz cuyo determinante es cero, Octave genera un aviso al 
tratar de calcular su inversa, aunque devuelve una matriz que resulta incorrecta: 


»C = [2 1; 
C = 

4 2] 

2 1 


4 2 


> det(C) 


ans = 


0 


>inv(C) 


warning: inverse: matrix singular to machine precisión, rcond = 0 

ans = 


4.00000 

2.00000 

0.50000 

0.00000 
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Ejercicios - 

S ► 17 Calcula los determinantes y, en su caso, las matrices inversas, de las siguientes matrices: 


'2 1' 


4 

-3 

2 ' 


0.1 

-5.0 

3.0 

8.7' 

4 2 

b) 

8 

5 

— 12 

9 

-5 

-2 

c) 

2.0 

2.7 

-1.6 

9.2 

4.5 

5.6 

8.9 

6.7 







0.2 

-4.5 

-8.0 

1.0 


H ► 18 En la siguiente tabla se muestra el coste por hora de cuatro procesos de manufactura 
distintos. También aparecen el número de horas de cada proceso necesario para fabricar tres 
productos diferentes. 


Proceso 

Coste/hora (€) 

Horas 

necesarias por 

unidad 

Producto 1 

Producto 2 

Producto 3 

Torneado 

10 

6 

5 

4 

Afilado 

12 

2 

3 

1 

Molienda 

14 

3 

2 

5 

Soldadura 

9 

4 

0 

3 


Escribe un programa Octave que, 

a) determine el coste unitario de cada uno de los productos, 

b) calcule el coste total de producción de 10, 5 y 7 unidades de los productos 1, 2 y 3, 
respectivamente. 

Deberás utilizar únicamente productos de vectores y matrices. 
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2.7. Aplicación 

A priori, la idea para resolver la aplicación planteada al principio del tema consiste 
en calcular con las fórmulas que se indican los valores de [A], [B] y [C] para distintos 
valores del tiempo. Ahora bien, ¿hasta que valor del tiempo debemos calcular las 
concentraciones?. 

Teniendo en cuenta que, según la reacción, el reactivo A va a desaparecer, podemos 
pensar que basta con esperar a que, por ejemplo, sólo quede el 1 % de A. Es muy fácil 
calcular a qué tiempo corresponde esta concentración. A partir de la ecuación para la 
concentración de [A], tenemos: 

[A] = 0.01 * [A]o — 0.01 * [A]o = [A]o * e ~ kut — —> t max = 

Con este valor de t max , resulta sencillo desarrollar el siguiente programa Octave 
para resolver el problema planteado. 


% Programa : cineticaABC.m 

7 0 Descripción: Este programa calcula y representa la variación de 
7o las concentraciones con el tiempo en una cinética 

7o del tipo A->B->C, con constantes kl y k2. 

7o Pedimos al usuario los valores de A0, kl y k2. 

A0 = input(’Introduzca el valor inicial de la especie A: ’); 

kl = input(’Introduzca el valor de la constante kl: ; 

k2 = input(’Introduzca el valor de la constante k2: J ) ; 

7o Calculamos el tiempo máximo y el vector de tiempos 

tmax = round(-log(0.01)/kl) ; 
t = linspace(0,tmax,100) ; 

7o Calculamos los valores de las concentraciones de A,B y C. 

A = A0*exp(-kl*t) ; 

B = kl*A0*( exp(-kl*t) - exp(-k2*t))/(k2-kl) ; 

C = A0*(1 + ( kl*exp(-k2*t) - k2*exp(-kl*t))/(k2-kl) ); 

7o Impresión de resultados 
dispC’t A B C') 
disp([V A ’ B ’ C’]) 

Ejercicios - 

19 Escribe este programa en un fichero de nombre cineticaABC.m y ejecútalo con los 
valores que se indican en la aplicación. 
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2.8. Ejercicios prácticos 

Es conveniente que pienses y realices los ejercicios que han aparecido a lo largo de la 
unidad marcados con el símbolo A antes de acudir a la sesión de prácticas correspon¬ 
diente. Deberás iniciar la sesión realizando los ejercicios marcados con el símbolo H. 
A continuación, deberás hacer el mayor número de los ejercicios siguientes. 

Ejercicios - 

► 20 Crea un vector x con los elementos x n = jjy para n = 1,2,3,..., 20 y comprueba 
que i B -¡0a medida que n aumenta. 

► 21 Dada la matriz A = [2 41; 672; 35 9] , escribe las órdenes necesarias para 

a) asignar la primera fila de A a un vector de nombre xl. 

b) asignar las 2 últimas filas de A a una matriz de nombre y. 

c) calcular la suma de cada una de las columnas de A. 

d) calcular la suma de cada una de las filas de A. 

► 22 Dados los vectores x = [1 4 8] e y = [2 1 5], y la matriz A = [3 16; 5 2 7], 
determina cuáles de las siguientes órdenes se ejecutarán correctamente en Octave y da su 
resultado. En caso de que la expresión dé un error, indica por qué. 

a) x + y 

b) x + A 

c) x’ + y’ 

d) [x ; y’] 

e) [x ; y] 

f) A - [x ; y] 

► 23 A partir de la matriz A = [2 797; 3156; 812 5], explica los resultados de 
las siguientes expresiones: 

a) A> 

b) A(:, [1 4]) 

c) A( [2 3] , [3 1]) 

d) A(1:3,:) 

e) [A; A(1 :2, : )] 

f) sum(A) 

g) sum(A’) 

h) prod(A) 

i) prod(A J ) 

j) sum(sum(A) ) 

k) mean(A) 

l) mean (A 1 ) 

m) mean(mean(A)) 
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► 24 Define las variables m = 10, n = 6, r = 3ys = 6,y construye las siguientes matrices: 

1. Matriz identidad de orden n. 

2. Matriz con elementos todos igual a —1 de orden s x r. 

3. Matriz cuadrada con elementos de la diagonal igual a 0 y el resto 1 de orden s. 

► 25 Define los siguientes vectores: 

vi = [1 2 3 • • • 10 ] v2 = [20 18 16 • • • 2] 

Analiza el comportamiento de la siguientes operaciones diciendo de qué tipo de operación se 
trata: vectorial, escalar, matricial, elemento a elemento, ... También indica las que producen 
error indicando por qué: 


a) 

vi + v2 

f) 

vl/v2 

k) 

sin(ir 

* ul). * v2 

b) 

vi * v2 

g) 

vl./v2 

1) 

(vV). 

* (v2') 

c) 

vi' * v2 

h) 

vi A 2 




d) 

vi. * v2 

i) 

vi. A 2 

m) 

( vi' * 

v2) A (-1) 

e) 

vi. A v2 

j) 

2 * (vi + v2) 

n) 

((«!') 

. * v2). A (-1) 


► 26 El factor de compresibilidad z mide la desviación del comportamiento de un gas res¬ 
pecto a la ecuación de estado de los gases ideales. Si el gas fuese ideal, z sería igual a la unidad 
bajo todas las condiciones, pero lo más habitual es que este factor se desvíe considerablemente 
de este valor. Además las desviaciones de la conducta ideal pueden hacer que z sea mayor o 
menor que la unidad. 

Una posible forma de calcular el factor de compresibilidad es utilizar la ecuación de estado 
del virial, según la cual: 


z = l + 


ai 

V 


a 2 03 

y2 + y3 + ■■■ 


Las cantidades ai, 02 , etc, se denominan el segundo, tercero,etc., coeficientes del virial, y 
dependen sólo de la temperatura y de las propiedades de las moléculas del gas. 

Escribe un Programa que pida al usuario los coeficientes del virial y el volumen de un 
gas y calcule el factor de compresibilidad z. 

Nota: Puedes escribir el programa sin utilizar vectores o con ellos. Puedes pedir al usuario 
los datos que consideres oportunos. Ten en cuenta que no asumimos a priori el número de 
coeficientes del virial. 


► 27 La media armónica es otra forma de calcular la media de un conjunto de N números. 
Esta media se calcula con la siguiente expresión: 

media armónica = -¡-,-=— 

— + — H-h — 

Xl X2 XN 

Escribe un programa Octave que lea un número arbitrario de valores positivos y calcule 
su media armónica. En caso de que aparezca algún número negativo o nulo, el programa 
terminará con un mensaje de error. Puedes elegir el método que quieras para pedir los datos 
al usuario. 
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